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More than 90% of the signal at an ocean-viewing, satellite sensor is due to 
the atmosphere (Gordon, 1987), so a 5% sensor-calibration error when viewing a 
target that contributes but 10% of the signal received at the sensor may result in a 
target-reflectance error of more than 50%. Since prelaunch calibration accuracies of 
5% are typical of space-sensor requirements, recalibration of the sensor using ground- 
based methods (e.g. Gordon, 1987; Carder et al. 1993) is required for low-signal 
targets. Known target reflectance or water-leaving radiance spectra and atmospheric- 
correction parameters are required. 

In this article we describe an atmospheric-correction method that uses cloud- 
shadowed pixels in combination with pixels in a neighboring region of similar optical 
properties to remove atmospheric effects from ocean scenes. These neighboring pixels 
can then be used as known reflectance targets for validation of the sensor calibration 
and atmospheric correction. The method uses the difference between water-leaving 
radiance values for these two regions. This allows nearly identical optical 
contributions to the two signals (e.g. path radiance and Fresnel-reflected skylight) to 
be removed, leaving mostly solar photons backscattered from beneath the sea to 
dominate the residual signal. Normalization by incident solar irradiance reaching the 
sea surface provides the remote-sensing reflectance of the ocean at the location of the 
neighbor region. 

Errors result from this method if corrections are not made for the following: 
1) reduction of the path radiance over the shadowed pixels resulting from partial 
shadowing of the overlying atmospheric column; and 2) adjacency effects resulting 
from solar-derived photons leaving the water from the bright, surrounding waters that 
are scattered into the field of view for the shadowed pixels. 

The initial application of the shadow-neighbor method was made assuming 
that the sensor calibration is correct. For illustrative purposes, assume that the solar 
zenith angle 0 Q is 45°, and that the angle from the pixel to the sensor 0 is 0°. A small 
cumulus cloud (about 1 km diameter) at 1.2 km altitude removes direct solar photons 
and shades a region. The water-leaving radiance from this shady region directed 
toward the sensor is designated L^ s . The wavelength dependence of terms in this 
article have been left out for brevity. This radiance results from skylight photons 
which have been scattered from beneath the ocean surface. 

Adjacent to the shadowed region is a neighboring patch of water with 
identical inherent optical properties to those of the shadowed region, differing only in 
the fact that direct, solar photons as well as skylight irradiate it. The water-leaving 
radiance from the neighbor region is designated L^. Since the cloud is quite small, it 
occludes a similar fraction of the sky for both shadow and neighbor patches, and 
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contributes a similar amount of diffusely scattered light to each. Furthermore, the 
atmospheric columns over each are considered to be identical. Therefore, both the 
path radiance contributed to the total signals and the diffuse attenuation t of the 
signals are nearly equal. Later, we will investigate the differences induced by the fact 
that a short column of the atmosphere over the shadow is shaded from solar photons 
by the cloud. 

To second order, the path radiance consists of a portion due to molecular or 
Rayleigh scattering, L r , one due to aerosol scattering, L a , and one due to multiple 
scattering, L ra . Let the radiance measured at the sensor when viewing an area in 
unshadowed water be given by (Wang and Gordon, 1992). 

0 Lin = L a + L r + + 1 Lwn 

where 

Lj n = measured radiance when viewing an unshadowed (neighbor) pixel 

L a = path radiance due to aerosol scattering 

L r = path radiance due to Rayleigh scattering 

L n = path radiance due to Rayleigh-aerosol interaction 

t = diffuse transmittance of the atmosphere for water-leaving radiance 

= water-leaving radiance of neighbor pixel 
L ws — water-leaving radiance of shadowed pixel. 

The apparent path transmittance of the water-leaving radiance from the 
shadowed pixel may not be equivalent to the term used in Eq. 1. Use of the "diffuse 
transmittance" is justified when viewing a large homogeneous area where the radiance 
scattered out of the viewing path is balanced by the radiance scattered into it from 
adjacent areas of the scene. In the case of the shadowed pixel, the adjacent areas of 
the scene are generally brighter, and so the apparent transmittance may be increased. 


If a terms express perturbations due to non-homogeneity in the scene 
illumination, the total radiance at the sensor when viewing a shadowed pixel may be 
written 


(2) L ls - L a - AL a + L r - aL f + L ra - AL ra + (t+ At) L ws . 

In general, the water-leaving radiance may be expressed as the sum of two 
parts: one part caused by backscattering of incident, diffuse skylight, and the other 
by backscattering of the direct solar beam by water and hydrosols. For the neighbor 
and shadowed pixels, respectively. 


(3) L wn sky^wn + 


sol^wn’ 


sky^w 


s + sol^ws* 


Under the assumptions that the shadowed pixel receives no direct sun and that 
the incident sky light at both pixels is the same, that is 


( 4 ) 

^ws sky^ws - skyL wn’ 

(5) L ln - L^ = aL 3 + AL r 4- a L ra + t so iL wn - At sky L wn . 

The first three terms on the right of Eq. 5 depend on the length of the 
shadowed portion of viewing path to the shadowed pixel. The height of the 
intersection of the viewing path and the upper edge of the cylinder of atmosphere 
shaded by the cloud can be determined from scene geometry. The layer of 
atmosphere below this intersection would normally be the source of less than 20% of 
the Raleigh scattering, so we will assume that the Raleigh-aerosol correction term, 
AL ra , must be negligible. Then, following Gordon et al. (1987), for the aerosol and 
Rayleigh corrections we may write 
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( 6 ) 

where 


aL x = {oj x t x FoP x (0 0 o> X)}t /4ir, x = a, r 


t x = optical thickness of shaded viewing path (vertical) 

< i) x — single-scattering albedo 

_ p^g-Cr-r’Vcostfo 

F 0 = extra-terrestrial solar irradiance 

P x (0 f 0 o ) = {P x (0_) + [p(0 ) + p(0 o )] P x (0 + )}/cos0 
cos 6 ± = ± cos6 0 cosd + sin0 o sin6cos(<j>-<j> o ) 

p(0) = Fresnel reflectance (0) 

P x (0) = scattering phase function 

t’ = diffuse transmittance from top of layer to top of atmosphere 

t’ = exp{-(r r /2 - r’/2 + r oz - r ra )/cos0}. 

The term involving At in equation (2) represents the apparent increase in 
diffuse transmittance of water-leaving radiance when viewing a shadowed pixel. The 
work of Tanre et al. (1979) indicates that this term is proportional to the difference in 
water-leaving radiances from the shadowed and neighbor pixels and is dependent on 
the geometry of the particular case. We will express At as 

(7) At = t * O (L^ “ — 1 * ^ (sol^wn ^sky^wn)* 

Preliminary Monte Carlo simulations indicate that a(550 nm) is of the order 
of 0. 1 for circular shadows of radius equal to 1 km for moderately turbid atmospheres 
(r a = 0.2). Thorough investigation of the dependence of a on wavelength and 
geometry is pending. Using equation (7) and equation (5) we obtain 

(8) solLwn = {Ltn ' Lis - AL a - AL r }/t(l-ff). 

Now define (Gordon et al. 1987) 

(9) w a(\) T a(\)P a (^ 0 ’\) 

£ (^i>\) ~ , 

a> a (tyT a (A j )P,(00 o >>j) 

then 

(10) AL a (Xj ) = € (X i ,\ j ) F o( X i) t frj) = S(Xj, Xj) 

a L a (Xj) f;(Xj) t (Xj) 

then 

(11) t(X i )[l-a(X i )] sol L wn (X i ) = L^fXi) - LJXj) - AL r (X ; ) - S(X i ,X j ) 

{L-t n (\j) - LJXj) - AL r (Xj) - t(Xj)[ l-cr(X j )] sol L wn (Xj)> 

If Xj is a wavelength such as 780 nm where the water-leaving radiance is 
essentially zero, then soiLw^Xj) = 0, and 

(12) aL 8 (780) = L, n (780) - L ts (780) - AL r (780) 

and 

(13) t(X i )(l-<r(X i )) sol L wn (X i ) = L. A) - MXi) - AL r (Xj) 

-S(X i ,780)[L ln (780) - L, s (780) - AL r (780)]. 

From scene geometry, the cloud height H and height of the shadowed 
viewing path can be calculated. Then, AL r (X) can be determined for all X by 
differencing Lowtran 7 results for the entire air column and the air column down to 
height H. The aerosol correction term, AL a (780), is then determined by Eq.12. 
Knowledge of wind speed allows estimation of €(Xj, 780) (Gregg and Carder, 1990), 
and thus S(\,780). The a(X t ) term involved in the transmissivity correction has been 
estimated by Monte Carlo simulations to be described later. Thus, Eq.13 allows 
^(L^Xj) to be calculated. 
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The results for an AV1RIS scene, 19 November, 1992 seaward of Key 
Biscayne, Florida are shown in Figure 1, where terms from Eq. 5 are normalized by 
the incoming solar irradiance (no skylight) to form remote-sensing reflectance 
spectra. Note that the spectrum resulting from the cloud-shadow method differs from 
that using conventional atmospheric correction methodology (see Carder et al., 1993) 
by only about 3.5%. This derived reflectance spectrum can in principle serve as a 
calibrated ground target to help with model parameterization for atmospheric 
correction of the entire scene, or in validation of the calibration of a sensor for scenes 
where conventional methods were unavailable. 
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Fig. 1: Remote-sensing reflectance spectra using conventional atmospheric-correction 
methods and the cloud-shadow method where path-shadow and transmittance 
perturbations are sequentially removed. 
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